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The invariant cross sections for direct photon production in hadron-hadron collisions are calculated 
for several initial energies (SPS, ISR, SppS, RHIC, Tevatron, LHC) including initial parton transverse 
momenta within the formalism of unintegrated parton distributions (UPDF) . Different approaches 
from the literature are compared and discussed. A special emphasis is put on the Kimber-Martin- 
Ryskin (KMR) distributions and their extension into the soft region. Sum rules for UPDFs are 
formulated and discussed in detail. We find a violation of naive number sum rules for the KMR 
UPDFs. An interesting interplay of perturbative (large k1) and nonperturbative (small fc?) regions of 
UPDFs in the production of both soft and hard photons is identified. The fct-factorization approach 
with the KMR UPDFs is inconsistent with the collinear approach at large transverse momenta of 
photons. Kwiecinski UPDFs provide very good description of all world data, especially at SPS and 
ISR energies. Off-shell effects are discussed and quantified. Predictions for the CERN LHC are 
given. Very forward/backward regions in rapidity at LHC energy are discussed and a possibility to 
test unintegrated gluon distributions (UGDF) is presented. 



PACS numbers: 

I. INTRODUCTION 

It was realized relatively early that the transverse mo- 
menta of initial (before a hard process) partons may 
play an important role in order to understand the distri- 
butions of produced direct photons, especially at small 
transverse momenta (see e.g.Q). The emitted photon 
may be produced directly in a hard process and/or from 
the fragmentation process. The latter process involves 
the parton-to-photon fragmentation functions which are 
not very well known. The isolation criterion used now 
routinly in the analysis of experimental data helps to re- 
duce the second component almost completely. 

The simplest way to include parton transverse mo- 
menta is via Gaussian smearing Q, |^ 1^. This phe- 
nomenological approach is not completely justified the- 
oretically. One should remember that there are differ- 
ent reasons for nonzero transverse momenta of incoming 
partons. First is purely nonperturbative, related to the 
Fermi motion of true hadron constituents. The trans- 



verse momenta related to the internal motion of hadron 
constituents are believed to be not too large, definitely 
smaller than 1 GeV. The second is of dynamical nature, 
related to QCD effects involved in the evolution of the 
parton cascades. The latter effect may be strongly de- 
pendent on longitudinal momentum fraction of the par- 
ton taking part in the hard (sub)process. 

The unintegrated parton distributions (UPDF) are the 
basic quantities that take into account explicitly the par- 
ton transverse momenta. The UPDFs have been studied 
recently in the context of different high-energy processes 
milillllllllSIIlIll- These works concentrated 
mainly on gluon degrees of fredom which play the dom- 
inant role in many processes at very high energies. At 
somewhat lower energies also quark and antiquark de- 
grees of freedom become equally important. Recently 
the approach which dynamically includes transverse mo- 
menta of not only gluons but also of quarks and anti- 
quarks was applied to direct-photon production 0|, 0| . 
In these calculations unintegrated parton distributions 
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proposed by Kimber-Martin-Ryskin were used. In 
this approach one assumes that the transverse momenta 
are generated only in the last step of the evolution ladder. 

Up to now there is no complete agreement how to in- 
clude evolution effects into the building blocks of the 
high-energy processes - the unintegrated parton distri- 
butions. In the present paper we shall discuss in detail 
a few approaches how to include transverse momenta of 
the incoming partons in order to calculate distributions 
of direct photons. 



two-dimensional function: 

foaussiki) = exp (-fcjV2o-o) /t^ • 

The UPDFs defined by Eq.|(2Jl and l|SJl are normahzed 
such that: 

j J^f — ^(x, kl^Ml) dkl = xpT'\x, h'f) . (4) 



II. UNINTEGRATED PARTON 
DISTRIBUTIONS 

In general, there are no simple relations between un- 
integrated and integrated parton distributions. Some of 
UPDFs in the literature are obtained based on famil- 
iar coUinear distributions, some are obtained by solv- 
ing evolution equations, some are just modelled or some 
are even parametrized. A brief review of unintegrated 
gluon distributions (UGDFs) that will be used also here 
can be found in Ref.jl^]. We shall not repeat all de- 
tails concerning those UGDFs here. We shall discuss in 
more details only approaches which treat unintegrated 
quark/antiquark distributions. 

In some of the approaches mentioned one imposes the 
following relation between the standard coUinear distri- 
butions and UPDFs: 



a(a;,yu^) 



fa{x. 



1 2 2 



'ir 



(1) 



where a = xq or a = xg. 

Since familiar coUinear distributions satisfy sum rules, 
one can define and test analogous sum rules for UPDFs. 
We shall discuss this issue in more detail in a separate 
section. 

Below we shall discuss in detail some of the approaches 
for UPDFs. Some other approaches are discussed e.g. in 
Ref.[l3. 



A. Gaussian smearing 

Due to its simplicity the Gaussian smearing of initial 
transverse momenta is a good reference point for other 
approaches. It allows to study phenomenologically the 
role of transverse momenta in several high-energy pro- 
cesses. We define simple unintegrated parton distribu- 
tions: 



\X, klfil) = Xpf\x, fil) ■ fGauss{kl) , (2) 



where ^^p) are standard coUinear (integrated) par- 

ton distribution (i = g, q, q) and fGauss(kf) is a Gaussian 



B. KMR distributions 



Kimber, Martin and Ryskin proposed a method to con- 
struct unintegrated parton distributions from the conven- 
tional DGLAP parton distributions 0|. Then 



fa{x,kf,^l'^) =Ta{kt,fl) 



as{kl\ 
2-K 



(5) 



where Paa'[,z) are splitting functions and a! {^^k\^ are 
parton densities, where a! = ^g or ^q. Angular-ordering 
constraint A = \ kt\) regulates the soft gluon sin- 

gularities. Recently Lipatov and Zotov [ll| used this 
method to calculate the direct photon spectra. Techni- 
cally they did not use the original KMR method. Instead 
they have written 



2tt 



dz 



P. 



{z)-q (-,ki) e(A -z) + P,,{z)-g ki) 



(6) 



fg{x,h,fi ) = Tg{k^,fJ. ) 



27r 



dz 



Y.P,,{z)-q (-,kA + P,,(z)^5 (^,kl) e(A - z) 



(7) 



In the following we shall call it LZ KMR prescription. 

The virtual corrections are resummed via Sudakov 
form factors: 



In Tq{kt,f?^ 



^ dp? a.(p^) 
2 p2 2tt 



dz Pqq{z) 



(8) 
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dPt a^(P?) 



2\ r 



2tt 



dz Pqg(z) 



dZ Z Pgg{z) 



(9) 



where Zmax = 1 - Zmm = m/(m + |Pt|)- 

The KMR method presented above can be used for 
transverse momenta fc^ > k^ Q. In the present paper we 
assume saturation of UPDFs for kf < k^Q. This is a bit 
arbitrary procedure. We shall discuss the consequences 
of the procedure for physical observables. 



C. Sum Rules for KMR UPDFs 

In order to gain more insight into the KMR distri- 
butions described shortly in the previous section in the 
following section we shall formulate and check some sum 
rules. 

Let us start from the valence number sum rules. We 
define the following integrals for up quarks: 
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^ dx 



dkf 



(10) 



and for down quarks: 



dki 



{x,kl^,^)~f^^^(x,kl^^^)] 



(11) 



The parameter fcj q is implicit for the KMR distribu- 
tions as discussed in the previous section. Naively one 
would expect: iV„^^, = 2 and Nd^^^ = 1. We shall check 
the dependence of these quantities on fj? and the freezing 
parameter fc^Q. The results are shown in Fig^ 

Somewhat surprisingly the results depend strongly on 
fj?' and the freezing parameter fc| q . The results are iden- 
tical for the standard KMR prescription and the one pro- 
posed by Lipatov and Zotov. In the integrals above the 
parameter /i^ occurs as an argument of the parton distri- 
butions as well as the upper limit of the internal integral. 
It seems interesting to allow for independent parameters 
in the two places. Therefore we define new quantities for 
up quarks: 



1 dx 



10 ^ 
'„2 ,2 



dkf 



{x,k^,iJ,Q) f^^^^{x,kf,fiQ)\ 



(12) 
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FIG. 1: The number sum rule as a function of fj, using 
standard KMR prescription with Sudakov form factor for up 
valence quarks (upper panel) and down valence quarks (lower 
panel) for several values of kt^. 



and for down quarks: 



NlJ^^',^,lkl,)^J^'^ 
x[fr'^{x,kl^^l)^ 



dkl 



ff''^^(x,kl^,l)] 



(13) 
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The results for several /Xq are shown in Fig[21 
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FIG. 2: The modified number sum rule as a function of fr^ 
using standard KMR prescription with Sudakov form factor 
for up valence quarks (upper panel) and down valence quarks 
(lower panel) for several values of po- 
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FIG. 3: The momentum sum rule as a function of /i^ using 
standard KMR prescription (upper panel) and Lipatov-Zotov 
prescription (lower panel) with Sudakov form factor. Here 
kio = 0.25 GeV^ 



Now a saturation of the sum rules for larger than These differences cancel as far as valence quarks are con- 
100 GeV^ can be observed sidered, which can be seen by inspection of Eq.© and/or 

Another interesting quantity is: Eq.®, Eq.(|7Il. 

xN^{^,^^,lkl,)^ fdx r [/™(.,fc^M^)] dkf 

Jo Jo 

(14) 

which can be interpreted as the contribution of parton of d. Kwiecinski unintegrated parton distributions 

a given type i to the momentum sum rule. In FigEl we 
show contributions for g, u, u,d,d,s — s as a function of 

the scale parameter . Kwiecinski has shown that the evolution equations for 

Again the defined above integrals are functions of the unintegrated parton distributions takes a particularly 

scale fi^. In this case there are huge differences for the simple form in the variable conjugated to the parton 

standard (left panel) and LZ (right panel) prescription. transverse momentum. In the impact-parameter space 
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the Kwiecinski equation takes the following simple form 



dfsix,b,fi^) ) 



2^ rl 



Pn 



, dz{Q{z - x) Jo{{l ~ z)fib) 



zPqqiz) + zPgq{z 



fs{x,b,i?) \ , 



la 



2^^2 / dzi^Qiz - x) Jo{{l - z)fib) 



Jiz)fs[^A^^')+Pggiz)fG{^,b,^i^) 
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fG{x,b,f? 



We have introduced here the short-hand notation 
Ins = fu — fu, fd ~ fd , 

fs = fu + fu + fd + fd + fs + fs ■ 



(15) 



(16) 



The unintegrated parton distributions in the impact fac- 
tor representation are related to the familiar coUinear 
distributions as follows 



(17) 



On the other hand, the transverse momentum dependent 
UPDFs are related to the integrated parton distributions 
as 



XPk[X,^J. 



dk't fk{x,kf,fi'^'' 



(18) 



The two possible representations, in the momentum 
space and in the impact parameter space, are interre- 
lated via Fourier-Bessel transform 



oo 



fkix,kt,fx ) = / dbbJo{ktb)fkix,b,^ ) , 
Jo 

poo 

fk{x,b,y?) = I dkt ktJo{ktb)fkix,kf,fi^) 
Jo 



(19) 



The index k above numerates either gluons (k=0), 
quarks (k> 0) or antiquarks (k< 0). While physically 
fk{x, kf , fj,^) should be positive, there is no obvious rea- 
son for such a limitation for fkix, 6, /i^). 



In the following we use leading-order parton distribu- 
tions from Ref. (231 as the initial condition for QCD evolu- 
tion. The set of integro-differential equations in b-space 
was solved by the method based on the discretisation 
made with the help of the Chebyshev polynomials (see 
jlTj'). Then the unintegrated parton distributions were 
put on a grid in x, b and and the grid was used in 
practical applications for Chebyshev interpolation. 

For the calculation of inclusive and coincidence cross 
section for the photon production (see next section) the 
parton distributions in momentum space are more use- 
ful. These calculation requires a time-consuming multi- 
dimensional integrations. Therefore an explicit calcula- 
tion of the momentum-space of the Kwiecinski UPDFs 
via Fourier transform for needed in the main calculation 
values of (xi^k^^) and (x2,A:|j) (see next section) is not 
possible. Therefore it becomes a neccessity to prepare 
auxiliary grids of the momentum-representation UPDFs 
before the actual calculation of the cross sections. These 
grids are then used via a two-dimensional interpolation 
in the spaces {xi, kf j) and {x2, fc| j) associated with each 
of the two incoming partons. 

The evolution of the parton cascade leads to a spread 
of the transverse momentum of the parton at the end of 
the cascade (the parton participating in a hard process). 
Let us define the following measure of the spread: 



^de,kUk{x,kiy^^) 



Jdkf fk(x,klfl^) 



(20) 



Above fk can be either gluon (fc = 0) or quark (fc > 0) or 
antiquark (k < 0) distribution. As an example in Fig^ 
we show the spread, obtained for different parton species, 
as a function of parton longitudinal momentum fraction. 

In this calculation the factorization scale was fixed at 
/i^ — 100 GeV^. The Kwiecinski evolution leads to in- 
creasing spread with decreasing longitudinal momentum 
fraction. The spread for different species of partons is 
quite different. In the region of small x the spread in k^ 
for gluons is bigger than a similar spread for sea and va- 
lence quarks. This is very different than a corresponding 
spread for Gaussian distributions which is usually taken 
to be independent of x and parton species. In addition, 
the spread of k^ for small values of x is considerably larger 
than the nonperturbative spread of the initial Gaussian 
distributions, taken here identical for all species and en- 
coded in the model parameter bo- 

In contrast to the KMR unintegrated valence quark 
distributions the Kwieciiiski valence quark distributions 
fulfill the number sum rules for up: 



1 dx 



X JQ 

Kwtec 



dkl 



X [f^^'^'ix, kl /i^) - fi-'''{x, kl^,')] 



(21) 
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E. Mixed distributions 



^ 10' 
cn 

o 10 

> 



1 



' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' 

Kwiecinski UPDFs 
/i' = 100 GeV^, bo = 1 GeV' 




J I I I I I I L. 



-4 -3 -2 -1 

log,o(x) 



FIG. 4: < fc? > as a function of x for different unintegrated 
parton distributions: solid - glue, dashed - Uvai, dotted - d^ai 
and dash-dotted - sea. 



The calculation with the Kwiecinski distributions 
discussed in the previous section shows that the spread 
in fcj for gluons can be much bigger than that for 
quarks/ antiquarks. On the other hand in the region of 
small X there are several unintegrated gluon distributions 
available in the Hterature. At high energy (small x) 
the contribution of qg^ gq subprocesses is larger than 
the contribution of qq, qq subprocesses. Therefore it 
seems reasonable to use the different UGDFs from the 
literature together with the Gaussian distributions for 
quarks and antiquarks as discribed above. Such an 
approach is especially justified at forward (y-y 3> 0) 
and backward photon rapidities (y^ ^ 0), where xi <C 

I, X2 ^ I {{kl,t)glue > {k2,t)q,q) OT Xl 1, 3:2 < 1 

{{kl,t)qq < {k2,t)glue)- 



III. UPDFS AND PHOTON PRODUCTION 



and for down quarks 



dkf 



122) 



The cross section for the production of a photon and 
an associated parton (jet) can be written as 



da{hih2 —>■ ^f, parton) 
(Ppi,t<Pp2,t 



= / dyidy- 



d'^ki^t (Pk2,t 



\Q'k'^{xiX2s)^ 



\M{i] ^ 7/c)|2 



d'^{ki,t + k2,t -Pi,t -p2.t)Mxi,kl Afj(x2,klt) , 



(23) 



where fci t and k2,t are transverse momenta of incoming 
partons. The longitudinal momentum fractions are cal- 
culated as 



Xi 



mit 



m2t 



{q,q,9), {q,q,9), {q,9,q), {9,q,q), etc. (see FigEJ. 

If one makes the following replacements 



(24) and 



fi{xi,klt) xiPiixi)S{klf) 



(26) 



X2 



mit 



m2t 

7^' 



(25) 



where mu and m2t are respective transverse masses. 
In the leading-order approximation: {i,j,k) = 



fj{x2,klt) X2Pjix2)S{klt) 



(27) 



then one recovers the standard coUincar formula (see 

The inclusive invariant cross section for direct photon 
production can be written 



FIG. 5: The diagrams included in our fet-factorization ap- 
proach with the notation of kinematical variables. 



7]- 'P2,t=ki^t+k2,t-pi,t 



(28) 



I 



and analogously the cross section for the associated par- 
ton (jet) can be written 



can be performed by introducing the following new aux- 
iliary variables: 



da{hih2 k) 
dy2d'^P2,t 



= [dy,^^ {...)]. , 

TT TT Ipi,t='£l,t+fe2, 



t-P2,t 

(29) 



Qt = + k2,t , 
qt = ki,t - k2.t ■ 



The integrand I{ki^t, k2.t',yi,Pi.t) defined in Ea. H28|l de- 
pends strongly on UPDFs used. We shaU return to this The jacobian of this transformation is 
interesting issue. 

Let us return to the coincidence cross section. The 
integration with the Dirac delta function in H23I) 



d{Qt,qt) 



1 1 

1 -1 



1 1 

1 -1 



2 • 2 



, , d'^ki t d^k2,t , N j2/ ^ 
dy\dy2 (...)(5^(...) 



(30) 



Then our initial cross section can be written as: 



(31) 



(32) 



da{h\h2 —>■ ^,parton) 1 



d'^Pi^td'^P2-, 



dyidy2 d^Qtd^qt (...) 5'^{Qt - Pi,t - P2,t) 



i J dyidy2 (fqt^ (••■) Iq.^p, 



7 / dyidy2 qtdqt dip (...) \g p 



^ / dyidy2 ^dqf dip (...) \q^^p^ 




r 



(33) 



Above Pt — pi,t + P2,t- Different representations of the bution of the sum of transverse momenta of the outgoing 
cross section are possible. If one is interested in the distri- 
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particles (parton and photon), then it is convenient to 
write 

cPpi,t d^P2,t = ^d^Ptd^Pt = ^dip+PtdPt dcp^ptdpt 



= - 2TTPtdPt dip^ptdpt . 



(34) 



If one is interested in studying a two-dimensional map 
Pi,t X P2.t then the differential volume element can be 
written 



d^pi.t d^p2.t = Pi,tdpi,t d(t>2 P2.tdp2,t ■ 
Then the two-dimensional map can be written as 



(35) 



da{pi^t,P2,t) 
dpi,tdp2.,t 



d(t)id(t)2Pi,tP2,t I dyidy2^qtdqtc 



u {■■■)■ 
(36) 

The integrals over 0i and (f)2 must be the most external 
ones. The integral above is formally a 6-dimensional one. 
It is convenient to make the following transformation of 
variables 

(01, 02) — > (0siim = 01 + 02, 0<ii/ = 01 - 02) , (37) 

Explicit formulae for the basic matrix elements are given 
in Appendix B. 







FIG. 6: Integrand under the ki^t and k2,t integration in the 
invariant cross section formula for proton-proton scattering at 
-^s = 63 GeV and y— 0, pt,-y = 5 GeV and Gaussian UPDFs 
((To = 1 GeV). Off-shell matrix elements for gluons are used, 
(a) qq, (b) qq, (c) gq, (d) qg. 



IV. INCLUSIVE PHOTON SPECTRA 

A. Integrands of the inclusive cross sections 

Before we go to the discussion of the dependence of the 
invariant cross sections on the values of rapidity and pho- 
ton transverse momentum let us consider the integrand 
J(fci_t, ^2,*; ?/i,Pi,t) (before integration over ki^t and k2.t) 
in Eq.^. 

In FigEl we show an example for ^/s = 63 GeV, y = 
and pt,-y = 5 GeV. In this calculation the unintegrated 
parton distributions based on GRV coUinear parton dis- 
tributions 22\ and Gaussian smearing (ctq = 1 GeV) in 
parton transverse momenta were used. This is a rather 
standard method to "improve" the coUinear approach. 
We do not need to mention that this, although having 
some physical motivation, is a rather ad hoc procedure. 
The two-dimensional distributions are peaked for small 
values of ki^t and k2^t- How fast the distribution decreases 
with ki^t and/or fc2,t depends on the value of the smear- 
ing parameter ctq- The larger ctq the slower the decrease. 

In Fig|3we show similar maps {^/s — 63 GeV, y — 0, 
pt .^ = 5 GeV) for the KMR UPDFs. Three local maxima 
can be seen in the figure. A first maximum occurs when 
both ki^t and k2^t are very small. This is caused by the 
structure of UPDFs themselves. The two other maxima 
occur when ki t = Pt and k2^t is small or fc2,t = Pt and 
ki^t is small. These are caused by the structure of matrix 
elements. The presence of long tails in kt in the KMR 




FIG. 7: Integrand under the ki^t and fc2,t integration in the 
invariant cross section formula for proton-proton scattering 
at ^ = 63 GeV and y= 0, pt,-y = 5 GeV and KMR UPDFs 
{kt,o = 0.25 GeV^). Off-shell matrix elements for gluons are 
used, (a) qq, (b) qq, (c) gq, (d) qg. 



distributions is a necessary condition to produce the sec- 



9 



ond and third maxima. When pt increases the second and 
third maxima move towards larger ki t and/or A:2,t- This 
clearly shows that the range of integration must depend 
on the value of photon transverse momentum. In Fig|S| 
we show some integrand (KMR UPDFs) but for larger 
energy W — 630 GeV and larger photon transverse mo- 
mentum pt^j = 50 GeV. 
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FIG. 9: Integrand under ki^t and ^2,* in the invariant cross 
section formula for proton-proton scattering at ^/s = 63 GeV 
and y= 0, pt,-, = 5 GeV and Kwiecinski UGDF (bo — 1 
GeV~^ , fj? = 100 GeV^). Ojf-shell matrix elements for glu- 
ons are used, (a) qq, (b) qq, (c) gq, (d) qg. 



B. Off-shell effects 



FIG. 8: Integrand under the k\^t and fc2,t integration in the 
invariant cross section formula for proton-antiproton scatter- 
ing at = 630 GeV and y= 0, pt,^ = 50 GeV and KMR 
UPDFs {klo = 0.25 GeV^). Off-shell matrix elements for 
gluons are used, (a) qq, (b) qq, (c) gq, (d) qg. 



Let us quantify the kinematical off-shell effect by defin- 
ing the following quantities: 



joff 



— shell 



TO f f — shell 



jon — shell 



(38) 



Fig|H| looks very much the same as Fig[7| if the trans- 
verse momenta of incoming partons are rescaled by the 
ratio of photon transverse momenta. 

In Figl^lwe show a similar map for the Kwiecinski dis- 
tributions. In this case the factorization scale is fixed for 
= 100 GeV2. 

All kinematical variables are exactly the same as in 
the previous cases. The integrand is rather similar to 
the one for the KMR UPDFs, except that the first max- 
imum at fci^t, fc2,t « is somewhat broader. We shall see 
consequences of the different integrands when discussing 
transverse momentum dependence of the photon inclu- 
sive cross sections. 



In Fig^Jwe present results for Rgq (left panels) and 
Rqg (right panels) for proton-proton scattering at W = 
63 GeV {y-y — 0, p-y^t = 5 GeV) and proton-antiproton 
scattering at W = 630 GeV {y-y = 0, p^,t = 50 GeV). 
When ki^t,k2,t the off-shell effects dissapear, i.e. 
the ratio becomes unity. The larger transverse momenta 
of gluons the larger the off-shell effect is. Therefore one 
may expect a related enhancement of the photon inclu- 
sive cross section when the UGDFs with large transverse 
momentum spread are used. 

In Fig^]we show the ratio of the inclusive cross sec- 
tions obtained with off-shell and on-shell matrix elements 
as a function of photon transverse momentum. In this 
calculation the Gaussian distributions were used with dif- 
ferent values of the ao parameter. The bigger ao the 
larger the enhancement due to the off-shell effects. 

In Fig ll2l we show similar enhancement for a few rep- 
resentative UPDFs discussed in section 2. The biggest 
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FIG. 10: Off-shell effects for the dominant mechanisms for 
pp scattering at ^/s = 63 GeV, y = 0, pt = 5 GeV (upper 
panels) and pp scattering at y/s = 630 GeV, y = 0, pt = 50 
GeV (lower panels). 



enhancement is obtained with the KMR and BFKL dis- 
tributions, i.e. those which have the biggest gluon trans- 
verse momentum spread. In general, the bigger photon 
transverse momentum, the smaller the enhancement. We 
conclude that at larger photon transverse momenta one 
can use standard on-shell matrix elements. 



C. Photon transverse momentum distributions 

Let us start the analysis from the lowest energies. In 
Fig ll3l we show inclusive invariant cross section as a func- 
tion of Feynman xf for several experimental values of 
photon transverse momenta as measured by the WA70 
collaboration. It is well known that the coUinear ap- 
proach (dotted line) fails to describe the low transverse 
momentum data by a sizeable factor of 4 or even more. 
Also the fct-factorization result with the KMR UPDFs 
(dashed line) underestimate the low-energy data. In con- 
trast, the Kwiecihski UPDFs (solid line) describe the 
WA70 collaboration data almost perfect. In order to il- 
lustrate the situation in Fig^]we show the ratio 



da 



Kwiec 



dy<Ppt 



da' 



coll 



dyd'^p, 



(39) 



as a function of Feynman xf for pt — 4.11 GeV and 
Pt = 5.70 GeV. This figure shows that the enhancement 
factor strongly depends on xf- Such an enhancement 
is required by the experimental data as can be seen by 
inspection of the previous figure. 
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FIG. 11: The ratio of inclusive cross sections with off-shell 
and on-shell matrix elements with Gaussian UPDFs and dif- 
ferent values of the parameter ctq as a function of photon 
transverse momentum. 



In Fig ll5l we show invariant cross section for direct 
photon production as a function of photon transverse mo- 
mentum for photon rapidity y — and ^/s — 63 GeV. 
The results obtained with the KMR UPDFs strongly de- 
pend on the value of the parameter fcj q. The larger the 
parameter fc^g, the smaller cross section. This means 
that even at large photon transverse momenta the non- 
perturbative effects (small fct's) play an important role. 
This can be better understood via inspection of the two- 
dimensional maps ki^t X ^2,* shown in Fig[71and is related 
to the second and third local maxima which give a sig- 
nificant contribution to the invariant cross section. 

In principle, one could try to find the parameter fc^g 
by confronting the theoretical results with experimental 
data. If the parameter is adjusted to larger transverse 
momenta there is a deficit at smaller transverse momenta 
compared to the ISR data . 
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FIG. 12: The ratio of inclusive cross sections with off-shell 
and on-shell matrix elements for different UPDFs as a func- 
tion of photon transverse momentum. 



In Fig llfil we compare our results with recent proton- 
proton RHIC data [2^. Here only low transverse- 
momenta of photons were measured. The results ob- 
tained with the Gaussian UPDFs strongly depend on the 
value of the cto parameter which is not surprising for the 
low transverse momenta. The Kwiecinski distributions 
give fairly good description of the PHENIX data. 

In contrast to the "low energy" data, there is no deficit 
for the KMR UPDFs at larger energies as can be seen 
from Figs. [13 dl and d The KMR UPDFs, however, 
strongly overestimate the experimental data at large pho- 
ton transverse momenta. This is especially visible for 
proton-antiproton collisions at W = 1.96 TeV when com- 
pared with recent Tevatron (run 2) data [33- FigstTHI 
and 1191 show that the unintegrated parton distribution 
approach with the KMR UPDFs is clearly inconsistent 
with the standard coUinear approach at large transverse 
momenta. This is caused by the presence of large-Zct tails 



^'yX ' ' 'Vs = 2i.'GeV' 
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Colllneor 




FIG. 13: Invariant cross section for direct photons for ^/s 
= 23 GeV as a function of Feynman xp for different bins 
of transverse momenta. In this calculation of shell ma- 
trix elements for subprocesses with gluons were used. The 
Kwiecinski UPDFs were calculated with the factorization 
scale fj,^ — 100 GeV^. The theoretical results are compared 
with the WA70 collaboration data H. 
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FIG. 14: The ratio of the invariant cross sections as a func- 
tion of If for ^/s — 23 GeV. The Kwieciiiski UPDFs were 
calculated with the factorization scale fj,^ — 100 GeV^. 
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FIG. 15: Invariant cross section for direct photons for ^/s = 
63 GeV and y = as a function of photon transverse momen- 
tum. In this calculation off-shell matrix element for gluons 
were used. The experimental data of the R806 collaboration 
are taken from Ref. |2!tI|. (a) Gaussian smearing (ao = 1, 2 
GeV) versus coUinear approach (b) standard KMR prescrip- 
tion, 



FIG. 16: Invariant cross section for direct photons for 
^/s = 200 GeV and y = Q. \n this calculation off-shell matrix 
element for gluons were used. The experimental data of the 
PHENIX collaboration are taken from Ref. |26|l . (a) standard 
KMR prescription (b) Gaussian smearing (ao = 1, 2 GeV) 
versus Kwieciiiski and coUinear approach 



(of the l/kt type) in the KMR UPDFs. It is not the case 
for the Gaussian and Kwiecinski UPDFs which seem to 
converge to the standard coUinear result at large photon 
transverse momenta. In this respect the latter UPDFs 
seems preferable. 

D. Direct photons at LHC 

Up to now we have confronted our results with the 
existing experimental data from SPS, ISR, RHIC, SppS 
and Tevatron. In a not too distant future one may expect 
experimental data from a newly constructed LHC. 

In Fig l20l we present transverse momentum dependence 
of the invariant cross section for the proton-proton col- 



lisions at W = 14 TeV for different values of photon ra- 
pidities. Different UPDFs lead to quite different results 
at the LHC energies. Therefore future measurements at 
LHC should give a chance to verify different UPDFs dis- 
cussed here as well as others to be constructed in the 
future. 

In Fig 12] we display average values of a;i and X2 of 
partons participating in a hard subprocess for two values 
of the photon transverse momentum pt = 10, 50 GeV. 
At large rapidities either xi 3> X2 or xi <^ X2- Then 
one expects the dominance of q{valence)g or gq{valence) 
hard processes (see Fia l^ . At such small values of x the 
evolution effects for UGDFs are expected to be very im- 
portant. In addition, one expects rather small transverse 
momenta of large-x valence quarks, much smaller than 
transverse momenta of the associated small-x gluons. 
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FIG. 17: Invariant cross section for direct photons for 
^/s — 630 GeV. In this calculation off-shell matrix element 
for gluons were used. The experimental data of the UA2 col- 
laboration were taken from Ref. 27] . (a) Gaussian smearing, 
(b) quark/antiquarks: Gaussian smearing (ao = 1 GeV), glu- 
ons: KL, BFKL and GBW, (c) standard KMR prescription. 
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FIG. 18: Invariant cross section for direct photons for 
^/s — 1800 GeV. In this calculation off-shell matrix element 
for gluons was used, (a) standard KMR prescription, (b) 
quarks/antiquarks: Gaussian smearing (ctq = 1 GeV), glu- 
ons: KL, BFKL and GBW. 



FIG. 19: Cross section for direct photons for ^/s = 1.96 TeV. 
In this calculation off-shell matrix element for gluons were 
used, (a) standard KMR prescription, (b) Gaussian smearing 
((To = 1,2 GeV) versus Kwiecinski UPDFs. 
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FIG. 21: Average values of xi (solid) and X2 (dashed) as a 
function of photon rapidity. The borders of the bands shown 
correspond to pt — 10 GeV (lower limit) ant pt — 50 GeV 
(upper limit). 
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FIG. 20: Predictions for LHC (^ = 14 TeV) for different 
photon rapidities y = 0, 2, 4. We compare results for differ- 
ent UPDFs: KMR(fc?,o = 1 GeV^), GRV+Gaussian smear- 
ing (o-o = 1 GeV) and BFKL, KL and GBW UGDFs with 
quarks/antiquarks smeared by a Gaussian distribution (cro — 
1 GeV). 




FIG. 22: The dominant mechansisms of photon production 
at forward (left panel) and backward (right panel) rapidities. 
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FIG. 23: The rapidity dependence of the invariant cross sec- 
tion for a few values of photon transverse momenta {pt = 10, 
20, 50 GeV). The contribution of gq (solid) and qg (dashed) 
subprocesses are shown separately for each value of the trans- 
verse momentum. 



Results obtained with different UGDFs differ signifi- 
cantly in the region of large rapidities. 

In the case of the Kwiecihski UPDFs at small trans- 
verse momenta {pt = 10, 20 GeV) only limited part of the 
full curve is shown. This is dictated by a purely techni- 
cal cut on the longitudinal momentum fraction x > 10~^ 
when constructing interpolation maps of the Kwiecihski 
UPDFs. As can be seen from Fig|^ large \y\ require 
very small xi or X2, sometimes smaller than 10~^. Fur- 
thermore the Kwiecihski distributions are not expected 
to be reliable at such small values of longitudinal momen- 
tum fractions. Therefore the range of application of the 
Kwiecihski UPDFs at LHC is limited in rapidity and pho- 
ton transverse momentum. The larger photon transverse 
momentum, the broder range of application in photon 
rapidity. 

In conclusion, the region of large rapidities {\y\ > 3), 
discussed in this section, seems an appropriate place to 
test models of UPDFs. It is not clear to us if any 
of the LHC detectors can register the large-energy for- 
ward/backward photons. The CASTOR detector associ- 
ated with the CMS detector is a potential option. 



V. CONCLUSIONS 

The inclusive cross section for prompt photon produc- 
tion has been calculated for different incident energies 
from SPS to LHC within the formalism of unintegrated 
parton distributions. Different models of UPDFs lead 
to rather different results. The Kwiecihski distributions 
provide the best description of experimental data in the 
broad range of incident energies. The existing exper- 
imental data test UPDFs down to x = 10~^, i.e. in 
the region of intermediate longitudinal momentum frac- 
tions adequate for application of the Kwiecihski equa- 
tions. Inclusion of the QCD evolution effects and espe- 
cially their effect on initial parton transverse momenta 
allowed to solve the long-standing problem of theoretical 
understanding the low energy and low transverse momen- 
tum data for direct photon production. 

As a buy-product we have analyzed momentum sum 
rule for different UPDFs. We have found that the KMR 
UPDFs violate naive number sum rules. The same dis- 
tributions lead to an interesting interplay of soft (small 
gluon fct's) and hard (large gluon kt's) regions of UPDFs. 
Even at large photon transverse momenta this interplay 
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causes a huge enhancement as compared to colhnear ap- 
proach, quite inconsistent with the experimental data at 
large photon transverse momenta. 

We have presented predictions for LHC based on sev- 
eral UPDFs with special emphasis on large rapidity re- 
gion. Here different UPDFs lead to quite different pre- 
dictions. Therefore we conclude that this region can be 
very useful to test different UPDFs. 
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VII. APPENDIX A 



which gives 



ki,tdki^tk2,tdk2^t = 



-J 10^iln(10)d2/il0^^1n(10)dy2 
i f 10^y'+y'hn\l0)dyidy2 . 



The moving with center-of-mass hadron-hadron energy 
maxima for the KMR distributions cause that the inte- 



Then the invariant cross section for the production of the 
photon (associated with a parton) can be writen as 



J 



o-j™(?7i,_Pi,t) 



d-qid'^pi,t 

- I d(t>id(t)2 ■ 10'°si°(''^i-'Hn(10)dlogio(fc? ,) • 10'°sio(fc2,t)in(lo)dlogio(fc2 J • drj: 
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-f,{xi,klt) fj{x2,klt) 



In the formula above fi{xi, k\ ^) and fj{x2, fc| J are un- cross section can be formally written as 
integrated parton distribution functions. The invariant 



J 



cTinv{m,Pi,t) = / dlogt^Q{klt)dlogiQ{klt) //og_zog(logio(A:J t),logio(A:| J) 



VIII. APPENDIX B 



The on-shell as well as off-shell matrix elements 



\M 

y^7-f P taken into account for the following sub- 
processes 

qq 73 (on — shell) , 

99 ~^ 19 {on — shell) , 

gq 75 {on — shell , o f f — shell) , 

19 ~^ 11 {on — shell ^ o f f — shell) . 



When neglecting parton masses the on-shell matrix ele- 
ments squared can be written as 
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\Maa- 


—f-'J^ 1 




\Mqq- 




= 7rQ;emVQ^l,sa2,s(167r) 






= TTttem V0^1,sa2,s(167r) 




-.7XP 


= 7raemVo^i,sa2,s(167r) 



9 y V i ii 
1\ / u s 



s u 
l\ ft s 

s i 



I 

Including finite mass efi^ects for quarks/antiquarks and written as: 
ofl-shell efi^ects for gluons |2fljl the matrix element can be 





-.7XP 




-.7XP 




-.7XP 



8(47r)2Q!g„^Q;i,sa2,s , , 

^o7? 2^27^^ ^^qq(M.t,k2.t) , 

8(47r)2Q!emV"l.s"2,s , , 

-777-^ iTTTT ^^qq[M,t,k2,t) , 

(47r)2Q!e„i^Q;i,sa2 
3(s-to2)2(^_^2)2^99 

em 

^^2\27I~Z272 



Fgq{ki,t, k2.t) 



\M,,^.X\^ = 3(,_^2)2(,':^2)2^-(fcM,^2.), 



r 



where for brevity we have introduced 



Fgg{ki,t,k2,t) = 6m^ - (3P + Sfi^ + 14i?i)m^ 

+ (P + ti^ + TtM^ + 7Pu^)ml 

- iu{P +y?) , 
Fqq{ki^t,k2,t) = &ml~{iu'^ + 2,P + lAui)m\ 

+ {u^ + P + 7uP + 7u^P)ml 

- ui{u^ + P) , 

Fgq{ki,t,k2^t) = 6m^~(2/cft + 2(su)fc^t + 3j2 + 3fi2 + 14su)m^ 

+ {2{su)kfi + Ssukfj. + + + Tsu^ + ls^u)m1 

- su{2k\ ^ + 2{su)klj^ + f + u^) , 
Fqgiki,t,k2,t) = 6m^^ - {2klt + 2{st)kl^ + Srs^ + 3P + Usi)m^g 



(2(s<)^_j + 8st/fc| t + .s-^ + + 7sP + 7sH) 
st(2A:2_t + 2(si)kl^ + f + P) 
I 



m 



q 



In the formula above only transverse momenta of the in- trix elements. Usually gluons generated via QCD effects 
going gluons are included explicitly when calculating ma- 
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have on average larger transverse momenta than quarks. 
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